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Cavity  Dynamics  for  Long  Rod  Penetration 

Minhyung  Lee  and  Stephan  Bless 


ABSTRACT 

Two  analytical  models  for  the  crater  size  due  to  long  rod  penetrations  are  presented.  The 
energy  method  is  the  first  which,  in  a  steady  state  penetration,  relates  the  kinetic  energy  loss 
of  a  penetrator  to  the  total  energy  deposited  in  the  target.  An  equation  of  cavity  dynamics  is 
set  up  by  assuming  a  radial  motion  in  the  target.  This  approach  appears  to  provide  an 
upper  bound  of  the  crater  size  due  to  the  neglect  of  energy  losses  in  the  penetrator. 

The  second  method  is  based  on  the  observations  that  two  mechanisms  are  involved  in 
cavity  growth  due  to  long  rod  penetrators;  flow  of  penetrator  erosion  products,  which 
exerts  radial  stress  on  the  target  and  opens  a  cavity,  and  axial  inertia  of  the  target  as  it  flows 
around  the  penetrator  nose.  The  analysis  includes  the  centrifugal  force  exerted  by  the 
penetrator,  radial  inertia  of  the  target,  and  the  strength  of  the  target.  Thus,  the  extent  of 
cavity  growth  due  to  penetrator  mushrooming,  which  cannot  be  predicted  by  other 
analyses,  can  be  estimated. 

The  analytical  models  are  compared  with  existing  formulas  and  empirical  equations  and 
found  to  be  more  accurate  than  other  analyses  in  the  literature.  Numerical  calculations 
using  AUTODYN-2D  are  also  conducted  and  these  confirm  the  phenomenological 
assumptions  in  the  models. 


I.  INTRODUCTION 

The  cavity  dynamics  and  the  resulting  crater  size  created  by  long  rod,  kinetic  energy  (KE), 
anti-armor  penetrators  are  examined  analytically  and  numerically. 

The  radial  motion  of  target  material  generated  by  hypervelocity  penetration  occurs  in  two 
stages  as  shown  in  numerical  simulations  [1].  In  the  first  stage,  eroded  penetrator  elements 
play  the  dominant  role  in  opening  a  cavity.  In  the  second  stage,  the  inertia  deposited  in  the 
target  is  responsible  for  the  further  cavity  expansion  until  the  strength  of  the  target  forces  it 
to  come  to  rest. 

An  earlier  analytical  study  of  the  crater  size  from  long  rod  penetrators  is  presented  by 
Szendrei  [2]  based  upon  a  momentum  argument.  By  accounting  for  the  first  stage  with  a 
proper  choice  of  initial  conditions,  Shinar  et  al.  [3]  set  up  an  equation  for  processes 
occurring  during  the  second  stage  and  determined  the  final  crater  size.  However,  there  are 
some  difficulties  associated  with  the  initial  conditions,  since  they  assume  that  the  inertia 
phase  starts  when  the  cavity  radius  is  equal  to  the  initial  penetrator  radius.  Based  on  an 
energy  balance  argument,  Tate  [4]  provided  an  approximate  formula  for  the  crater 
diameter.  Bjerke  et  al.  [5]  presented  an  empirical  equation  for  the  channel  diameter  versus 
impact  velocity  by  fitting  a  polynomial  to  data  for  tungsten-nickel-iron  sintered  penetrators 
ranging  from  90  to  94%  tungsten  impacting  thick  and  normal  incidence  steel  targets. 
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In  the  present  work,  we  propose  two  analytical  models  for  the  crater  size.  The  analytical 
models  are  compared  with  existing  formulas  and  empirical  equations.  An  energy  method 
states  that  the  kinetic  energy  of  a  penetrator  is  converted  to  the  plastic  work  for  crater 
production  and  kinetic  energy  deposited  in  the  target.  This  approach  can  provide  the  upper 
bound  of  the  cavity  size  due  to  the  neglect  of  energy  absorbed  during  deformation  of  the 
penetrator. 

The  second  model  uses  the  method  of  Miller  [6]  to  determine  the  amount  of  the  first  stage 
cavity  expansion  due  to  eroded  penetrator  element.  Because  of  the  neglect  of  the  target 
strength  in  his  model,  the  cavity  radius  is  not  bounded.  In  the  present  work,  including  the 
strength  of  the  target  allows  us  to  explore  the  ultimate  cavity  radius  achieved  during  the 
first  stage  cavity  expansion,  due  solely  to  a  mushrooming  effect.  We  then  use  the  energy 
method  in  order  to  estimate  the  amount  of  the  second  stage  cavity  expansion  due  to  inertia. 

In  both  models,  an  axisymmetric  coordinate  system  moving  at  the  constant  penetration 
velocity  U  is  used  as  reference.  Numerical  simulations  using  AUTODYN-2D  are  also 
presented  in  comparison  with  the  analytical  predictions. 


II.  ENERGY  PRINCIPLES 


The  mathematical  formulation  of  the  cavity  dynamics  for  penetration  of  long  rod 
penetrators  is  described  in  this  section  within  the  framework  of  an  energy  balance  principle. 
The  basic  assumptions  required  in  the  analyses  are: 


—  Target  material  is  considered  as  incompressible  and  inviscid. 

—  The  motion  of  the  cavity  wall  is  purely  radial. 

—  The  internal  energy  change  due  to  thermal  effects  and  strain  energy  in  the 
target  are  negligible. 

—  Energy  absorbed  by  the  projectile  is  negligible. 

—  Steady  state  is  assumed. 

—  The  penetrator  material  ultimately  comes  to  rest. 

A.  Kinetic  Energy  Loss  of  a  Penetrator  Per  Unit  Target  Length 


Consider  a  long  rod  penetrator  of  density  p^,  diameter  Dp,  and  length  L  impacting  an 
infinite  target  of  density  p^  at  normal  incidence.  For  a  penetrator  with  an  impact  velocity 
V,  the  kinetic  energy  loss  of  the  penetrator  per  unit  target  length  is  given  by. 


F  = 


f) 


(1) 


where  m  =  p  (^r  ^).  Here  m  is  the  mass,  P  the  final  penetration  depth  and  r  the 


2 


radius  of  the  penetrator.  Using  the  hydrodynamic  relation  L/P=  (V  -  U)AJ,  where  U  is  the 
penetration  velocity,  Eq.  (1)  becomes. 


E 


kP  ~ 


(V:U)v2 

u 


(2) 


With  this  result,  an  energy  balance  equation  will  be  set  up  for  the  purpose  of  determining 
an  equation  of  cavity  dynamics. 


B.  Equation  of  Cavity  Dynamics 


It  is  assumed  that  the  kinetic  energy  loss  of  the  penetrator  per  unit  target  length  equals  the 
total  energy  deposited  in  the  target  (kinetic  energy  plus  plastic  work  done  to  open  a  cavity). 
If  the  radid  velocity  is  u(r),  then  following  [7],  the  source  strength  is  related  to  the  radial 
distance  r. 


u  =  ^/r.  (3) 

The  kinetic  energy  per  unit  target  length  E^,  in  the  radial  field  within  a  finite  radius  Q  is 
then. 


E 


kt 


n 


j(7tp,  u^t)dT  =  Kp,g^  N  , 


(4) 


where,  N  =  ln(Q/a)is  a  dimensionless  geometric  parameter  and  a  is  a  cavity  radius.  If 
we  replace  Q  by  na  and  n  is  assumed  to  vary  in  the  range  5-15,  N  can  be  treated  as  a 
constant.  In  the  calculation  N  equals  2.3.  The  amount  of  plastic  work  done  to  open  the 
cavity  per  unit  target  length  is  given  by, 

a 

Ep,  =  J(2;rr)R,  dr  =  ;za'R,  ,  (5) 

0 


where  R,  is  the  target  resistance  for  radial  cavity  expansion.  Now  the  conservation  of  the 
energy  per  unit  target  length  can  be  expressed  by  the  relation. 

E,p  =  Kp^g'^  N  +  R, .  (6) 


By  introducing  constants. 


^kp 


;rR, 


we  can  obtain  an  expression  for  the  source  strength  given  by, 

^  =  QVp'-a'  . 


(7) 


(8) 
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In  Eq.  (3),  the  source  strength  is  also  ad.  In  order  to  obtain  the  equation  of  cavity 
dynamics,  Eqs.  (3)  and  (8)  are  combined  to  give, 


ad  =  Q^|?^-a^  .  (9) 

Integration  of  Eq.  (9)  with  a  boundary  condition  u  (t  =  0)  =  0  yields, 

a(t)  =  VP'-(P-Qt)^  .  (10) 


Equation  (10)  determines  the  time  evolution  of  the  cavity  radius.  The  cavity  has  a 
maximum  at  time,  t„j„  =  P/Q.  Note  that  Eq.  (10)  is  valid  only  for  0  <  t  <  t^^ .  At  t^^^, 
the  velocity  field  vanishes  throughout  the  flow  field  so  the  kinetic  energy  in  the  target 
becomes  zero.  The  maximum  cavity  radius  is  then  equal  to. 


By  substituting  Eq.  (2),  we  have , 


a 


c 


R, 


(V-U) 

u 


(12) 


A  similar  formula  to  Eq.  (12)  can  be  determined  from  momentum  principles  and  given  by 
(Appendix  A), 


a 


C 


I  ^  2p,(y~^ 

i  R,  ^  R, 


(13) 


where  Yp  is  the  dynamic  yield  strength  of  the  penetrator.  Tate's  formula  [4]  obtained  from 
an  energy  balance  argument  is. 


a 


C 


2p,(y-\jf 

R. 


(14) 


If  U=V/2  and  2pp(V-U)^  /Ri  >  1,  this  equation  coincides  with  Eq.  (12).  Equation  (13) 
has  the  same  form  as  Eq.  (14)  but  with  1  replacing  Yp  /R;.  Although  Eqs.  (12),  (13)  and 

(14)  are  similar  in  form,  they  are  derived  from  different  arguments.  Equations  (12)  and 
(14)  are  based  upon  transformation  of  penetrator  kinetic  energy  to  work  done  during  cavity 
expansion  and  Eq.  (13)  is  based  upon  momentum  exchange  between  the  penetrator  and 
target.  These  equations  will  be  compared  quantitatively  in  the  comparison  section. 
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C.  Determination  of  the  Penetration  Velocity  (U) 


For  steady  state  penetration  the  modified  Bernoulli  equation  applies, 

ip,(V-U)^+Y,  =  ip,U’  +  R,  .  (15) 


Then  U  is  given  by  [4], 


V-«^V^+  2(l-«')(R,-Yp)/p, 


(16) 


where  a  =  .  The  theoretical  R,  value  determined  by  the  classical  spherical  cavity 

expansion  model  as  presented  by  Satapathy  and  Bless  [8]  is  used  in  the  present  work: 


R. 


(Y  \ 
Jgj 


+ 


(17) 


which  gives  for  RHA  steel  the  value  Rj  =5.0  GPa . 


In  Figure  1,  the  present  model  prediction  for  the  evolution  of  the  cavity  radius  is  compared 
with  experimental  results  for  copper  jets  obtained  by  de  Rosser  and  Merendino  [9]  and  the 
analytical  results  given  by  Shinar  et  al.  [3].  The  impact  velocity  was  7.6  km/s  and  rp  was 

taken  to  be  1.75  mm.  The  yield  stress  of  their  steel  target  was  taken  to  be  0.8  GPa.  It  can 
be  seen  that  the  agreement  between  the  present  model  and  experiments  is  excellent  for  the 
evolution  period. 

Cavity  Radius  vs.  Time 


Figure  1.  Evolution  of  the  cavity  formed  in  RHA  steel. 
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III.  TWO  STAGE  CAVITY  EXPANSION  MODEL  (TSCEM) 


We  propose  that  in  the  first  stage,  the  eroded  penetrator  elements  exert  pressure  on  the 
target  and  open  a  cavity.  In  the  second  stage,  the  inertia  imparted  to  the  target  is 
responsible  for  the  further  expansion  of  the  cavity,  especially  during  hypervelocity 
penetration.  It  is  hypothesized  that  the  final  cavity  radius  is  equal  to  the  sum  of  the  cavity 
radii  produced  by  each  stage  acting  independently,  regardless  of  the  order  of  application  of 
loads.  In  other  words,  we  use  the  principle  of  superposition  for  the  total  work  done  by 
each  stage.  The  analysis  includes  the  centrifugal  force  exerted  by  the  penetrator,  radid 
inertia  of  the  target,  and  the  strength  of  the  target.  The  basic  assumptions  required  in  the 
analysis  are: 


—  Target  material  is  considered  as  incompressible,  steady  inviscid  flow  relative  to  the 
stagnation  region. 

—  The  pressure  profile  exerted  by  the  target  is  assumed  to  be  jp^U^sin^  p  as 
provided  in  Ref.  [6],  where  P  is  the  angle  between  the  tangent  direction  of  the 
centerline  and  the  axis  of  symmetry.  Although  the  Newtonian  pressure, 
jpJJ^sirP  P,  allows  a  simple  closed-form  solution,  it  over-predicts  the  pressure 
profile  [6].  This  pressure  is  acting  normal  to  the  flow  at  radius  r,  where  r  is  the 
radius  of  the  centerline  of  an  eroded  projectile  element. 

—  The  stagnation  region  over  which  the  pressure  equals  7 p,U^ exists  for  \  <P<k. 
The  domain  for  the  outside  of  the  stagnation  region  isO</J<;r/2. 

—  For  the  inertia  calculation,  we  make  an  analogy  with  the  case  of  a  rigid  body 
penetrator.  This  is  due  to  the  fact  that  a  projectile  head  moving  at  constant  velocity 
U  induces  a  kinetic  energy  transfer  to  the  target.  For  simplicity,  the  energy 
method  described  in  the  previous  model  is  then  adopted. 

A.  First  Stage  Cavity  Expansion 


With  these  assumptions,  the  first  stage  cavity  expansion  can  be  derived  from  Miller's 
analysis  [6]  with  presence  of  the  strength  of  the  target.  The  target  strength  will  eventually 
halt  the  cavity  growth.  The  geometry  of  the  first  stage  cavity  expansion  is  shown  in 
Figure  2. 


From  mass  conservation  in  a  coordinate  system  located  at  the  penetration  front,  the 
thickness  of  the  eroded  projectile  element  is  given  by. 


xiP)  = 


2r(i3)  ■ 


(18) 


Since  the  pressure  force  is  acting  normal  to  the  flow  at  radius  r,  the  centrifugal  forces  in  the 
eroded  element  for  the  outside  of  the  stagnation  region  must  be  balanced  by  the  resistance 
which  is  equal  to  the  pressure  force  plus  the  strength  of  the  target. 
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Figure  2.  Geometry  of  the  first  stage  cavity  expansion  model. 


R,  +  sin’;8  =  T(ft  p,  ,  (19) 

where  R(fi)  is  the  scaled  local  radius  of  curvature  of  the  center  line  curve.  in  this 
equation  is  approximately  an  average  of  the  spherical  and  cylindrical  cavity  expansion 
stresses.  From  geometric  relations. 


R()8)  =  - 


1  dr 
sinjS  d/3  ’ 


(20) 


Combining  Eqs.  (18),  (19)  and  (20)  gives  the  differential  equation  describing  the  trajectory 
of  an  eroded  projectile  element  given  by, 

R,  +  sin^P  =  (V-U)^  sinj3(^)-'  .  (21) 

2  2  r  dp 

By  introducing  two  constants. 


5P,(v-U)“ 


T  = 


R. 


Eq.  (21)  becomes. 


2rdr  = 


1  +  Tsin^^ 


(22) 


(23) 


In  order  to  determine  the  trajectory  of  the  eroded  element  as  a  function  of  P,  this  equation 
is  to  be  integrated  with  an  initial  condition,  r(P  =  7t  1 2)  =  which  yields. 


7 


r(/?)  = 


(24) 


il  +  Tsin')3 


By  integrating  Eq.  (24)  numerically  for  the  entire  domain  of  0  <  P  <  7r/2,  we  can 
determine  the  final  crater  size  (a^,).  Miller  [6]  made  a  hypothesis  that  there  is  a  stagnation 

region  in  which  the  pressure  is  equal  to  iPtU^  ,  such  that  r,  /tp  >  1.  Note  that  if  the  target 

pressure  is  approximated  to  be  i  sin'jS,  there  is  a  simple  closed-form  solution  given 

by, 


Ai  =  fp. 


v^py 


+ 


Vt  +  t^ 


ln[  1  +  2T  +  2Vt  +  T^ 


(25) 


In  order  to  find  the  initial  condition  rj ,  by  the  same  reasoning  used  above,  Eq.  (19)  can  be 
written  inside  of  the  stagnation  region  as, 


R,  + 


T()S)  Pp 


(V-U)^ 

R(i3) 


(26) 


Again  by  using  the  radius  of  curvature  and  the  thickness  of  the  elements,  this  equation  is  to 
be  integrated  with  an  initial  condition  r  (^  =  -jdl)  =  tp  /  2,  which  yields. 


r,  =  rpj--h 


2S 
1  +  T 


(27) 


B.  Second  Stage  Cavity  Expansion 

In  addition  to  first  stage  cavity  expansion  due  to  centrifugal  forces  exerted  by  the  eroded 
penetrator  elements,  the  inertia  imparted  to  the  target  is  responsible  for  fiirther  cavity 
expansion.  This  inertia  is  created  by  flow  in  the  target  around  a  penetrator  head  moving  at 
constant  velocity  U.  In  order  to  estimate  the  additional  cavity  expansion  due  to  target 
inertia,  we  consider  target  response  to  an  "equivalent  penetrator"-  a  rigid  body  of  the  same 
shape  penetrator  head  and  the  same  penetration  velocity.  We  can  then  use  the  energy 
method  described  previously. 

Consider  a  rigid  penetrator  moving  through  the  target  with  Poncelet  resistance.  Then  the 
force  (energy  per  unit  target  length)  may  be  written  as  follows, 

F  =  A  +  BU\  (28) 
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where. 


(29) 


A  =  (»i„")R„  B  =  l/2(m,.,")C,A  . 

Here  is  the  drag  coefficient.  In  these  equations,  A  is  the  target  resistance  force 
associated  with  plastic  work  and  B  is  the  drag  force  which  comes  about  due  to  the 
convective  inertia  effects.  For  the  second  stage,  this  inertia  is  responsible  for  the  cavity 
expansion  from  a^,  to  a^.  Hence, 

=  j2)lrR,dr=  .  (30) 

^cl 

The  solution  to  the  final  crater  size  is  then  given  by. 


=  a 


cl  ' 


1  + 


R. 


(31) 


The  drag  coefficient  is  assumed  to  be  0.5  for  a  hemispherical  nose.  The  extent  of  the 

compensation  due  to  inertia  increases  with  impact  velocity  because  the  term  in  Eq.  (3 1) 
becomes  more  important  in  the  high  impact  velocity  limit. 


IV.  NUMERICAL  SIMULATIONS 

Numerical  simulations  of  tungsten  long  rod  projectiles  (L/D=10)  into  infinite  RHA  steel 
targets  were  conducted  using  an  AUTODYN-2D  finite  difference  code.  The  objective  of 
the  simulations  was  to  confirm  the  phenomenological  assumptions  made  in  the  analytical 
models,  and  to  compare  numerical  results  to  the  analytical  results.  We  use  a  constant  shear 
modulus  G  and  a  von-Mises  yield  surface  with  a  yield  strength  Y  as  shown  in  Table  I. 


Table  L 

Shear  Stress 

Parameters 

Steel 

Tungsten 

G  (GPa) 

120 

140 

Y  (GPa) 

1.2 

2 

AUTODYN  is  an  Eulerian  wave  propagation  code  allowing  for  a  maximum  grid  of  60,000 
cells.  The  Mie-Gruniesen  equation  of  state  (E.O.S.)  is  used  for  both  materials.  The 
parameters  are  the  same  as  in  Ref.  [10]  and  listed  in  Table  n. 
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In  all  runs  a  constant  subgrid  with  eight  ceils  across  the  radius  of  the  projectile  and 
continuing  for  three  diameters  away  from  the  projectile  is  used.  The  penetration  direction 
has  a  constant  cell  size  to  make  the  subgrid  a  1 : 1  aspect  ratio. 

Table  11.  E.O.S.  Parameters 


Steel 

Tungsten  Alloy 

p.(g/cc) 

7.86 

17.4 

C  (km/s) 

3.57 

4.03 

S 

1.92 

1.26 

To 

1.7 

1.7 

P..  (GPa) _ 

-2. 

-2. 

V.  COMPARISON  OF  ANALYTICAL,  NUMERICAL,  AND 
EXPERIMENTAL  RESULTS 

The  crater  diameters  normalized  by  the  projectile  diameter  obtained  from  the  analytical  and 
experimental  results  versus  impact  velocities  are  shown  in  Figure  3.  The  results  are 
compared  with  the  numerical  predictions  obtained  by  using  AUTODYN-2D  with  eight  cells 
across  the  radius.  The  numerical  simulations  slightly  underestimate  the  crater  size,  which 
is  also  found  in  the  work  of  Bjerke,  et  al.  [5]. 

As  shown  in  Figure  3,  the  results  obtained  from  the  first  stage  cavity  expansion  analysis  is 
rather  good  until  V>2  km/s.  This  is  due  to  the  fact  that  the  inertia  effects  become  important 
in  this  range.  However,  the  analysis  can  provide  the  extent  of  the  mushrooming  versus 
impact  velocity.  By  including  inertia  in  the  two  stage  cavity  expansion  model,  the 
andytical  predictions  compare  well  to  the  data. 

It  should  be  emphasized  that  the  model  based  on  energy  conservation  provides  an  upper 
bound  for  the  final  crater  size  because  energy  losses  due  to  penetrator  deformation  are 
ignored.  That  is,  it  is  assumed  that  aU  of  the  available  kinetic  energy  of  the  projectile  is 
used  to  open  the  cavity.  It  can  be  shown  that  Eqs.  (12)  and  (14)  match  Tate's  formula 
quite  well  in  the  high  velocity  range  but  the  Tate’s  formula  appears  to  over-predict  crater 
size  in  the  low  impact  velocity  range.  The  agreement  between  the  energy  and  momentum 
models  is  quite  good.  The  equations  for  the  crater  diameter  obtained  from  different 
methods  are  shown  in  Appendix  B. 


VI.  CONCLUSIONS 

Two  analytical  models  for  the  crater  diameter  formed  by  long  rod  penetrations  are  presented 
and  compared  with  existing  formulas  and  empirical  equations. 

An  advantage  of  the  energy  principle  is  that  it  can  provide  an  upper  bound  for  the  crater 
size.  It  can  also  provide  the  evolution  of  the  cavity  radius.  The  formula  from  the 
momentum  principle  appears  to  provide  good  agreement  with  experiments  despite  the 
crudeness  of  the  model. 

An  estimation  of  the  amount  of  mushrooming  with  respect  to  impact  velocity  is  a  valuable 
contribution  of  this  analysis.  Maximum  penetration  for  a  given  kinetic  energy  should  occur 
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with  a  reduction  of  the  mushrooming,  reading  to  a  smaller  resistance  force,  as  discussed  by 
Magness  and  Farrand  [11]. 


-  Bjerke  (1992) 

-  Shinar  et  al.  (1995) 

■ 

AUTODYN 

- A — 

-FSCE 

— • — 

-  Energy 

- o  - 

-  Momentum 

—  ^  - 

-  Tate  (1986) 

- A  -  - 

- TSCEM 

L I ^ 1 ^ \ 1 I t I ! I ^ ^ 1 1 ^ I ^ \ j 

1  2  3  4  5 

V  (km/s) 


Figure  3.  Ratio  of  cavity  diameter  to  rod  diameter  versus  impact  velocity.  Tungsten 
alloy  striking  RHA  steel  target.  FSCE  represents  the  results  from  the  first  stage 
cavity  expansion.  AUTODYN  represents  results  obtained  with  the  AUTODYN 
wavecode.  Models  are  summarized  in  Appendix  B. 
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APPENDIX  A 


MOMENTUM  PRINCIPLES 


Consider  long  rod  penetration  into  a  homogeneous  target.  A  control  volume  attached  to  the 
penetration  interface  of  two  bodies  is  shown  in  Figure  A.  The  control  volume  is  bounded 
on  the  left  by  the  penetration  interface,  on  the  right  by  the  Hugoniot  elastic-plastic  interface, 
and  on  the  top  and  bottom  by  the  material  interfaces.  The  penetration  interface  is  moving 
with  velocity  U.  Mass  enters  the  control  volume  through  the  elastic  plastic  interface.  To 
determine  the  final  cavity  size  approximately,  the  conservation  of  momentum  equation  for 
the  control  volume  is. 


2p,(fD=)(V-U)“+  Y,(|d=)=  jR,(2ffl-)dr  =  R,(^D,')  ,  (Al) 


where  is  the  crater  diameter  and  is  the  radial  cavity  expansion  stress.  The 
conservation  of  mass  equation  for  the  control  volume  is  introduced  into  this  equation  by 
equating  the  incoming  mass  to  the  outgoing  mass  in  a  steady  state  condition.  The  solution 
for  the  final  crater  diameter  is  then  given  by. 


D  j 


R.  R. 


(A2) 


If  Yp  /R,  is  replaced  by  one,  this  expression  coincides  with  Tate's  formula  (Eq.  (14)). 


V-U 
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APPENDIX  B 


EQUATION  OF  CRATER  SIZE 

Energy  principle: 

Momentum  principle: 


Tate  [4]: 


TSCEM: 


Szendrei  [2]: 


Shinar  et  al.  [3]: 


Bjerke  et  al.  [5]: 


a,  =  rp(l.  1524  + 0.3388V +  0. 1286V')  . 


(12) 

(13) 

(14) 

(31) 

(Bl) 

(B2) 

(B3) 

(B4) 
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